Coherent Structures in Dense Granular Flows 
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We present a theoretical derivation of a rheology for dense granular flow, based on the process of 
inelastic collapse of neighboring particles. This collapse creates regions of correlated motion, which 
control the viscous behavior of the flow on a large scale. The result is a rheology that obeys the 
scaling form observed experimentally by Pouliquen and by the "G.D.R. Midi" group author. We 
identify the nature of the constraints imposed by inelastic collapse on the grain-scale motions in the 
flow; finally, we show using an energy cascade argument that the inelastic collapse need not proceed 
to the final endpoint in order for the correlations we have identified to build up. 

PACS numbers: 45.70.-n, 45.70.Mg, 45.70.Ht 
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Dense granular flows are ubiquitous in nature and in 
technology [1]. Notwithstanding, robust, microscopically 
well-founded descriptions of the behavior of these flows 
across flow geometries have proven elusive [2]. Several 
recent developments have hinted at the outline of such a 
description. In 1999, Pouliquen published a seminal pa- 
per performing a scaling analysis of "chute flows" on an 
inclined plane, showing that the rheology of these flows 
was controlled by a single scaling length [3]. These re- 
sults were reproduced in numerical studies [4] , and a po- 
tential explanation of this scaling length as arising from 
correlations in granular motion was advanced by the au- 
thors of this Letter [5, 6]. This scaling length has also 
been seen experimentally [7] . Most recently, a group au- 
thor "G.D.R. Midi" has published a comprehensive com- 
parison of numerical and experimental results on dense 
granular flows in a number of geometries, including both 
confined systems and systems with free surfaces [8] . They 
argue that for a system made of particles of diameter d, 
of density p, flowing at a shear rate 7, at a pressure P, 
the rheology is controlled by the local scaling variable 
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which figures as a generalization of the Pouliquen result. 

In this Letter we extend our earlier work on chute flows 
to account for Eq. (1) as a scaling variable for dense gran- 
ular flows. We present a microscopic picture of the origin 
of this quantity in the dynamics of the inelastic collapse of 
the particles in a flow, and we show how to apply these 
ideas to generate a systematic description of two stan- 
dard flows: pure shear at constant pressure, and chute 
flow with a free surface. 

It is well known that inelastic collapse is a key phe- 
nomenon in many types of granular systems [9]. Con- 
sider two particles brought together by a restoring force 
F at a normal velocity vq. If the particles have a resti- 
tution coefficient e, and they are perfectly rigid so that 
the collisions are instantaneous, we expect that they will 



collapse onto one another in a time 
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Real granular particles are not perfectly rigid, and 
may also have velocity-dependent coefficients of restitu- 
tion [10]. In practice, inelastic collapse will be controlled 
by a balance between a time scale r c that measures the 
duration of collisions and a time scale Tb that measures 
the intervals of ballistic motion between collisions. For 
perfectly rigid particles r c — -> 0, but of course Tb — > 
as well in the late stages of inelastic collapse. Thus a 
thorough analysis of inelastic collapse in a flow requires 
careful specification of the elastic and dissipative mecha- 
nisms of particle collision, as well as a clear sense of the 
flow conditions to which the particle is subject. 

We are first going to finesse this complexity by simply 
assuming that the time scale T nnc to reach the regime 
in which r c 3> Tb is short compared to the inverse shear 
time 7 _1 . Unfortunately, if one takes the limit of infinite 
rigidity of particles while keeping the coefficient of resti- 
tution and shear rate constant, one will never reach this 
regime. For this reason, at the end of this Letter we will 
consider the case in which Tb 3> t c always, and show that 
the results for the flow rheology are similar. 

Moving from the case of two particles to the case of a 
dense granular system at a pressure P, it is natural to 
generalize Eq. (2) to 
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where the local pressure P provides the restoring force, 
and where the unknown function /(e) accounts for the 
effect of multiple collisions with other particles in the 
course of an inelastic collapse event. The varying fre- 
quency of collisions needed to maintain the pressure on 
the grain may cause additional dependence of /(e) on the 
local packing fraction, which we have suppressed in this 
simplified analysis. 

For a granular system in motion with a local shear 
rate of 7, the time scale 7 -1 sets the time over which the 
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neighbors of a particle are changed by the overall shear 
flow. Provided 7T„ nc <C 1, we expect that clusters will 
form, the motion of whose particles are correlated with 
one another by this process of inelastic collapse. The for- 
mation of these clusters will be driven by the diffusive 
process of neighboring particles collapsing onto one an- 
other; each such collapse event takes a time of T nnc and 
introduces correlations into the flow on a scale of the 
interparticle distance, or the particle diameter d in the 
collapsed state. Thus the length scale of such correlations 
Ie will be determined by 

lE = Jfl, (4, 

y ' nnc J 

where d 2 /T nnc figures as a diffusion constant for a process 
occuring over a time scale of 7 -1 . We thus obtain 

V /( e ) p 7 

Note that there is no explicit dependence here on the 
particle diameter d. 

The precise nature of the correlations inside a clus- 
ter is somewhat subtle. When particles collapse onto 
one another, the relative normal velocity of the particles 
is driven to zero. In addition, the influence of friction 
will couple the rotational motion of the particles to their 
transverse velocity; immediately, in the case of infinite 
friction, or after a finite time, in the case of finite friction. 
Consider a set of N particles i in three dimensions with 
velocites v% and rotational velocities u)i (the latter the 
axial vector corresponding to the rotation about an axis 
cD). We further suppose the existence of pairs of particles 
(ij) that have inelastically collapsed onto one another, 
and for which there is no relative motion of the surface 
points in contact (corresponding to frictional locking of 
the particles.) Suppose that the vector connecting the 
collapse pair (ij) is Aw^jy Then the requirement that 
the relative motion of the surface points in contact be 
zero is 

Vi - Vj = Av {ij) =^(u)i + u)j) x Aw {ij) . (6) 

Taking the derivative with respect to time of this con- 
straint yields a constraint for the accelerations of the 
particles Si and angular accelerations = j^uji 

Aa {t]) = ^(wi+Wj) x Av^j) + ^(Ti + Tj) x Aw {ij) , (7) 

where, fortunately, complications owing to the non- 
Abelian nature of the rotation group do not appear 
to this order in d/dt. We call states constrained by 
Eqs. (6,7) "gear" states, for the obvious reason that the 
particles are rolling over one another like gears. 



These equations significantly constrain both the mo- 
tion of the particles and the forces between the particles. 
If the average coordination number of the collapsed par- 
ticles is z, then Eq. (6) gives three constraints per contact 
or 3z/2 per particle. Since each particle has 6 degrees of 
freedom without contacts, this means that the effective 
number of degrees of freedom per particle Ne/N is 

We thus see that the average coordination z < 4 in order 
for the collapsed state to be mobile at all, and Eq. (8) 
gives the effective number of degrees of freedom for the 
collapsed state [11]. 

The forces are even more completely determined. The 
total number of contact forces is exactly the same as 
the number of constraints on the accelerations given 
by Eq. (7), with the result that all of these forces are 
uniquely determined by the contact network and the ve- 
locities and angular velocities of the particles. 

Of course, the constraint equations on the velocities 
and accelerations do not hold indefinitely. At a certain 
point in the flow of the collapsed region, one of two break- 
downs in these equations must occur (see Figure 1). The 
first of these is the collision of two particles, which will 
initiate a new episode of inelastic collapse, at the end of 
which the contact network {(ij)} and the velocities will 
be altered by the various impulses created in the network 
of collapsed particles by this collision. The second type 
of failure is associated with the forces between particles. 
These are not arbitrary- the normal forces Af between 
granular particles may not be negative, and the tangen- 
tial forces T for particles with finite coefficients of friction 
k must obey a Coulomb constraint T < kj\f. As a pack- 
ing moves, we expect all forces determined by Eq. (7) to 
develop smoothly for some finite length of time, until one 
of these constraints is violated. At this time, a contact 
will fail and the contact network will change discontinu- 
ously. 

Having elucidated the nature of the correlations in the 
collapsed region, we return to the role of these regions 
in controlling the overall rheology of a granular flow. 
Our fundamental assumption is that the collapse scale 
Ie plays a role analogous to a Prandtl mixing length in 
hydrodynamics. Consider first the consitutive equation 
relating the shear component of the stress tensor a to the 
local shear rate 7 

a = fij, (9) 

defining a viscosity /i. We can write this viscosity in 
terms of a length scale 1^, 

p = pili, (10) 

effectively defining 1^. So far these manipulations are 
constrained only by dimensional analysis. We now intro- 
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stress differences develop [13]. This yields 




FIG. 1: Two means of disrupting the evolution of a collapsed 
coherent state: a) formation of a new contact at C due to 
inter-particle collisions-the particles are rolling on contacts A 
and B, and b) failure of an old contact at A under tension, or 
due to exceedance of a Coulomb criterion on the tangential 
force T < kN for normal force iV at a contact. 



duce the physical assumption that we can write 
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where this relation is written in the form of an asymptotic 
expansion in d/lE, and a, b are constants. This relation 
is the heart of our approach. Note that the sign of the 
first finite size correction bd/ls is constrained by the fact 
that conventional Bagnold scaling would imply that when 
Ie — > 0, Ifj, — > d [2]. (Although once we reach the limit 



d more conventional kinetic theory approaches may 



be superior to our approach.) 

We can now solve for the full shear dependence of p. 
From Eqs. (5,10), we see that 
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which together with Eq. (9) defines the rheology (ex- 
cepting more subtle effects such as normal stress dif- 
ferences) This is a Bingham-type rheology, with a yield 
stress added to a dynamic stress linear in the shear rate, 
although the pressure-dependences are novel compared 
to typical Bingham rheologies [12]. The alert reader will 
note the scaling dependence on the parameter / defined 
by the G.D.R. Midi group (Eq. (1)). 

Now let us consider the application of this result to 
some standard flow geometries. Consider first simple 
shear. In this case we fix the normal force per unit area 
F/v exerted into a sample by a boundary, as well as the 
tangential force per unit area Ft exerted at that bound- 
ary. Let us presume that we can write the pressure P as 
P = (1 + K)Fn/3, with K ^ 2 in the case where normal 



Ft 
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y<(l + K)F N /3p 



(14) 

showing a linear increase of the shear rate with tangential 
stress beyond a yield stress proportional to the normal 
force. Note that for a cell of size L, we must impose 
d < Ie < L in order for the inelastic collapse approach to 
make sense; indeed, we expect the flow not to be mobile 
if Ie > L, because collapse onto the wall will become 
possible on time scales shorter than 7 -1 . 

Turning to the case of flow down an incline at an angle 
6, for which a = pgsm9, and P = |(1 + K)pg cos 6*, we 
obtain 



tan 8 



3/(e) 



id 



+ K)g cos 6/3 



(15) 

It was shown already in Ref. [5] that this type of formula 
accounts well for the phenomenology of chute flows, as 
expounded in Ref. [3]. 

Our assumption that particles are correlated through 
collapse events of finite duration is somewhat naive-in 
practice the generation of "thermal" energy in the state 
through the impulses generated by the creation of new 
contacts (see Figure 1) will tend to disrupt contacts. It 
is therefore instructive to consider in more detail the en- 
ergy balance in a collapsing packing. In this way we can 
actually generalize our results to cases where inelastic 
collapse need not proceed to the end state of contacts of 
arbitrarily long duration. 

Our approach is to write formulae for the dissipation 
of energy on a succession of scales, and require that they 
be equal. Let us first consider the work done on a macro 
scale. From continuum theory, the rate at which work is 
done by the flow is 



V 



M = 07. 



(16) 



Returning to Figure 1, we see that in the collapsed state 
(which is our highest-level description of this process), 
the only dissipative process is that indicated in Figure la, 
in which a particle collides with another particle. Since 
the particles are part of extended objects of scale i M , the 
effective momentum difference in this collision is actually 
Ap oc mjl^, so that the rate of dissipation on a time scale 
where collapse events appear instantaneous is 



v c = pYil 



(17) 



Analogously to our discussion above, we can regard this 
as a definition of l^, but with the interpretation that this 
result corresponds to the dissipation associated with a 
particle collapsing onto an extended object of size ~ Z M . 

Now we resolve the flow on finer time and length scales, 
on which the particles are not fully collapsed onto one an- 
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other. We suppose that the particles are actually sepa- 
rated by small distances S -C d, with ballistic flight times 
across these distances (between collisions) of r b and col- 
lision durations of r c . If t c 3> Tt,, then in practice the 
packing is collapsed. This is the case that we examined 
above, which led to the fundamental rhcology given by 
Eq. (13). If ti, 3> t c , then the packing is somewhat sim- 
ilar to a granular gas; the difference between our ap- 
proach and normal kinetic theory approaches to motion 
is that in considering motion on scales S relative to the 
collapsed state, we effectively consider motion relative to 
the "gear" state defined by Eqs. (6-7). Normal kinetic 
theory approaches to flow consider thermal motion as 
being defined relative to a much simpler laminar flow in 
determining the "thermal" component of velocity [14]. 

If there is a constant coefficient of restitution of e, the 
loss of energy through inelastic "thermal" collisions is 



V 1 



Pi 1 - £ )Z3' 



(18) 



which is similar to Eq. (12), with the exception that the 
ratio I^/Ie appears to the second power in the former 
case, and to the first power here. Applying Eq. (11) 
connecting L to Ie, we see that this difference actually 
amounts at first order in d/ls only to a factor of two in 
the shear-rate dependent term in the effective viscosity, 
as well as changes to the overall constants multiplying the 
formula. Thus, while the precise nature of the inelastic 
collapse, and in particular the relative magnitudes of n 
and r c , do affect details of the rheology, they do not affect 
the qualitative result. 
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where the typical mean-free path S <C d in a dense gran- 
ular system. Clearly we require 



v M = v c = v I . 



(19) 



For a non-zero coefficient of friction, the transverse 
translational and rotational motion of the particles will 
lock as determined by Eq. (6) after a small number of col- 
lisions; we will get an approximately correct result if we 
assume that the particles lock after one collision, which 
will indeed be the case for an infinite coefficient of fric- 
tion. The size of a correlated region Ie is then determined 
by reprising the diffusive argument leading to Eq. (4), 



n t" 



(20) 



Again, we expect Eq. (11) connecting to Ie to hold. 
We can write the pressure from ordinary kinetic theory 
arguments as 



Poc 



1 mS 



or 



P_-d8 
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(21) 



(22) 



with B a constant, which completes our system of equa- 
tions [15]. Note that while the pressure is constrained by 
the order of magnitude expression Eq. (22), the average 
force between particles on time scales long compared to 
Tb but short compared to j^ 1 is still fixed by the gear 
equations Eqs. (6,7). 

After some manipulation, it is possible to reduce 
Eqs. (16-20,22) to 



°_ = M7 = V(l - e)^ 
P ~ P " ' Bl E 



(23) 
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